Spontaneous Lorentz Breaking and Massive Gravity 
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We study a theory where the presence of an extra spin-two field coupled to gravity gives rise to 
a phase with spontaneously broken Lorentz symmetry. In this phase gravity is massive, and the 
Weak Equivalence Principle is respected. The newtonian potentials are in general modified, but we 
identify an non-perturbative symmetry that protects them. The gravitational waves sector has a 
rich phenomenology: sources emit a combination of massless and massive gravitons that propagate 
with distinct velocities and also oscillate. Since their velocities differ from the speed of light, the 
time of flight difference between gravitons and photons from a common source could be measured. 

PACS numbers: 04.50.+h 



Since the early days of General Relativity (GR) there 
has been attempts to modify gravity at long distances. 
The problem is not just a theoretical challenge, but 
it could bring to important phenomenological conse- 
quences. The only ghost free massive deformation of 
linearized GR in the Minkowski background was given 
in [1] by Pauli and Fierz (PF). The peculiarity of grav- 
ity became apparent when it was realized that in the 
zero mass limit the PF theory was discontinuous [2] 
due the presence of an extra polarization state. The 
mystery of massive gravity deepened when it was re- 
alized that the propagation of five degrees of freedom 
is spoiled by interactions and a sixth ghost-like state 
turns on [3]. The problem was reexamined in the 
framework of effective field theory realizing that the 
reason behind the misbehavior of PF massive gravity 
was strong coupling [4]. It remains an open question 
whether a well behaved non-linear extension of the PF 
theory exists (see for instance [6]) but Lorentz invari- 
ance clearly plays a crucial restrictive role. If one gives 
up Lorentz invariance a trouble-free massive deforma- 
tion of gravity where the strong coupling and disconti- 
nuity issues are disentangled can be found [7, 8]. 

In this letter we describe a mechanism of sponta- 
neous lorentz breaking that provides mass to the gravi- 
ton in a consistent way. 

Let us consider a theory with two dynamical metrics 
92^v [15] each interacting with its own matter, 
a so called bigravity theory. The action contains two 
Einstein-Hilbcrt (EH) terms, and a mixed term [5]: 
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where £1.2 are the corresponding matter lagrangians. 

The mixed term V contains only non-derivative cou- 
plings between two metrics, therefore it can only be 
function of the tensor Xff — 9i ,ia g 2a ^[Q)- The cosmo- 
logical terms can be included in V, e.g. Va 1 = Aiq -1 / 4 , 
with q = detX — {g 2 j ' gi) 1 ^ ■ The equations of motion 
(EoM) are 
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where V" is the derivative of V with respect to Xjf. 
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The indices of the two equations are raised/lowered 
with the corresponding metrics. 

The action is invariant under a generic infinitesimal 
diffeomorphism (diff) generated by 

Sg a ^ = Q^Qaav + dvCdana + g*d a 9anv , (3) 

with a — 1,2. Notice that in the absence of V the sys- 
tem has a larger gauge symmetry: one can use different 
£f and £2 f° r 9i an d 92- The "diagonal" diff invariance 
is encoded in a set of generalized Bianchi identities. 

The two metrics can be diagonalized simultaneously, 
but in general their eigenvalues will not be propor- 
tional, and thus local Lorentz invariance will be broken. 
For vacuum solutions, we will assume that rotational 
invariance is preserved and that the two metrics have 
the same signature. 

The vacuum. The EoM (2) always admit constant 
curvature solutions for both g\ and g 2 , and in addition 
curvatures are proportional [10]. For simplicity here 
we focus only on the flat limit of those solutions, for 
which one fine tuning on V is necessary, as in standard 
GR for setting the cosmological term to zero. In this 
biflat case, the EoM are simply 



V = 0, 



V 



0. 



(4) 



where the bar stands for the background values. 

For rotationally invariant backgrounds, these are 
three independent equations. One of the equations cor- 
responds to the mentioned fine-tuning for biflat back- 
grounds. Then assuming that we live in sector 1 we 
can set 

gipy =diag(-l, 1,1,1) 

giixv = diag(-c 2 , 1, 1, 1) . 

Thus, the two remaining equations determine the con- 
stants c and u> for any given V. Physically c is the speed 
of light in sector 2, while u> parametrizes the relative 
conformal factor. However, a solution with c = 1 is al- 
ways present, since in this case two equations coincide 
(and determine to). 

Summarizing, we have two branches of solutions: 
Lorentz Invariant (LI) for c = 1, and Lorentz Breaking 
(LB) for 1 [16]. 

The LB branch is of particular interest since it nat- 
urally allows for consistent massive deformations of 
gravity. 
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Linearized analysis. Let us consider perturbations 
around the background (5), denned as g\ uu — g\^ v + 



In the LB phase, 1, this implies that Mi—0 and 
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The diffeomor- 



phisms (3) act at lowest order as 5h\ nu = 2g la ^ p d l/ ^ a 
and 8h2^u = 2u~ 2 g~2a(pd v )£, a . Since the background 
preserves rotations, we decompose both perturbations 
according to their spin content: 
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with <%u aj = 9^ = 9jXaij = SijXaij = 0. We have 
thus 1+1 gauge invariant tensors, 2+2 vectors, 4+4 
scalars. Also the 4 diagonal diffs can be split into one 
vector £ji and 2 scalars £°, <9j£\ As a result, we ex- 
pect that 2 tensors, 3 vectors and 6 scalars are physi- 
cal states determined by the EoM, while the remaining 
vector and 2 scalars can be gauged away. 

At quadratic level, in general the Lagrangian has the 



form L = Lhi 



kin 



+ L S 



where L^in contains 



derivative terms (in space and time) emerging from the 
EH actions, L mass comes the the quadratic expansion 
of the mixed term, while L source describes the gravita- 
tional coupling with matter. In field space, each field 
is a 2-component column vector, and it is convenient 
to define the following 2x2 matrices: C — diag(l,c), 
M 2 = M 1 2 diag(l, k) where n = M$/Mfu) 2 c. Then the 
kinetic term has the compact form 

Lkin = \x%M 2 (C 2 A - d 2 ) Xij - \wtM 2 AWt + 



2$*M 2 At-t*M 2 (C 2 A-3<9 2 )t , (7) 



where W l — u l — d t Si, $ = ^ — 2d t v + d 2 (j. Notice that 
Xij, Wi, $, t arc all gauge invariant fields making the 
invariance of L^n manifest. When the matter energy- 
momentum tensors are conserved also the source la- 
grangian is expressed in terms of gauge invariant fields: 
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The crucial term, L mass , in the flat limit stems from 
the expansion of V. For the rotational invariant bifiat 
background (5), for which V = V' = 0, it reads 

L mass = -2(3!5 2 ) 1/4 Tr (X\ V" X\) 

= \ (M a Q % Q h b m + 2MfK t h b M - MfK 3 h h l3 

+ Mf^ i h%-2Mth a m ht i ), (9) 

where X\ — Xg^h^ — g^hY X is the fluctuation of X. 
The second line defines the Lorentz breaking masses 
M [7] that can be computed explicitly once a V and 
a consistent background are given. Their pattern is 
drastically different in the two branches. 

Gauge invariance (3) gives crucial constraints: 

M>,l( c V)=0, Mi.2,3,4 Q =0, 

(1 c 2 )M 4 = 0. (10) 
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In the LI phase instead c = 1 and two conditions in (10) 
coincide, thus allowing a non- vanishing Mi. Also, all 
the A^'s are proportional to the same projector (11). 
Moreover, in this branch the mass term reduces to the 
generalized PF form: L mass = rf- a r\ v ^h t iiv Ah <x ^-\-h i 'Bh 
(A and B being 2x2 matrices). This is equivalent to 
M = A + B, M lt2 = -A, M 3A = B. We remark 
that the limit c — > 1 is discontinuous. 

Let us analyze in detail the linearized theory, sepa- 
rately in the two branches. 

Lorentz invariant phase. Since here all the mass 
matrices are proportional to the same projector with 
a null eigenvalue, the linearized theory is diagonalized 
by the mass eigenstates 
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where h cl ' 2 are the canonically normalized fields and 
the mixing angle is sin# = (1 + k) -1 / 2 . The state h c+ 
is massless, while h c ~ has a generalized PF mass term. 

The massive sector h c ~ will be plagued by the known 
problems related to the Pauli-Fierz case. In the gen- 
eral case, with A ^ —B, it contains ghosts [17]. If 
we remove the ghosts setting A = —B, we face the 
vDVZ discontinuity problem [2]. In this theory how- 
ever, since the original gravitons have different Planck 
scales, the discontinuity in our sector can be controlled 
by the mixing angle. The present bounds [13] would 
require M 2 ~ 30Mi. In addition the Newton constant 
becomes distance dependent: it changes by a factor of 
1 + 4 tan 2 9/3 from short to large distance, the critical 
distance being the inverse graviton mass. The matter 
of type 2, if it exists, will interact with ours in distance 
dependent way [14]. 

Finally, the sixth mode [3] of h c ~ will manifest once 
interactions are considered, and in general will propa- 
gate in both gravities, leading to strong coupling [4]. 

Lorentz breaking phase. We give an overall analysis 
of the static and propagating modes in the LB branch, 
and refer to table I for the full summary of all phases; 
details will be given in [10]. 

a. Tensors Each tensor has two independent com- 
ponents. They propagate according to the EoM: 



(k 2 -k 2 ) 








M 2 



K (k 2 -c 2 k 2 ) M 2 



Xij=0, (13) 



describing two gravitons with different "speeds of 
light" . The gravitons are mixed by the non-diagonal 
mass matrix M2, and thus oscillate and have a non- 
linear dispersion relation. Since M 2 has a zero eigen- 
value, expanding the dispersion relation in powers of 
k 2 one finds at low energy a massless graviton that 
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TABLE I: Summary of phases in the biflat limit. Propagating states are labeled as 1 2p | , |2 m | for spin 2 massless, massive. 
GR± refers to standard GR potentials for the ± combinations of <£>, r. GRdisc refers to GR potentials plus discontinuity, 
GR/i refers to GR potentials plus linear term (14). (*) Here fi 2 — A2(3As — A2)/(A3 — A2)M 2 . ( 2 ) A flat direction of the 
V potential, 5(i/>_,t_) oc (3A4,Ao), or Weyl_ for 3A4 = — Ao: these guarantee the vanishing of the extra linear potential. 
( 3 ) The source Tqq of the undetermined field should vanish; the orthogonal combination <$>d has Newtonian potentials. 



travels with speed v 2 = (1 + c 2 k)/(1 + k), and a mas- 
sive one, of mass m 2 = (1 + k~ 1 )\ 2 /M 2 . In the high 
energy limit two states propagate with different speeds 
and also oscillate. 

Interestingly, when c > 1, the second graviton prop- 
agates faster than "our" light, though this will not 
lead to causality violations. Indeed, in the coordinates 
where the metrics are (5), the Cauchy problem is glob- 
ally well posed in terms of the preferred time t. In prin- 
ciple this scenario could be tested by observing the time 
of flight difference between gravitational waves and op- 
tical signals, or frame dependence of the gravitational 
waves propagation that would provide the evidence for 
a preferred frame. 

b. Vectors Thanks to Mi = 0, the vector states 
do not propagate. The EoM however determine three 
static potentials: the two Wu, Wu that couple to the 
respective T i's as in standard GR, without mixing, 
and the combination — (Su — S 2 j), that is zero. 

c. Scalars Also scalar modes do not propagate, 
again because of Mi = 0; they however mediate the 
static potentials. By defining X 2 = A 2 + A (A 2 — A3), 

A 2 = 3Ai+A (A 2 -3A 3 ) and Tqo = (T 100 -T 200 /c^ 2 k), 
the solutions of the EoM are 
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(15) 



The two remaining gauge invariant fields, ip- = (c 2 xpi — 
^2) and o- = ((Ti — cr 2 ) are determined in terms of 
t_ = (ti — r 2 ). Since they have no source, the explicit 
expression is omitted. From <I>i in (14), Mi is identified 
as the Planck mass, M\ = M 2 p = (87rG) _1 . 



Since in this phase only gravitons propagate, strong 
coupling problems [4] arc absent. 

Linear term. In (14) we recognize the standard GR 
potentials plus, in $1,2, a term proportional to 1/A 2 , 
that represents a linearly growing potential at large 
distances [11]. For generic Ao,2, 3 ,4, the scale at which 
this force sets out is proportional to fi, (15). Since this 
linear growth signals the breakdown of perturbation 
theory at large distances from sources, it is interest- 
ing to study the conditions under which fi vanishes. 
One, cheap, possibility is to have massless gravitons, 
A 2 = 0; the other possibility is that the numerator 
vanishes, A 2 = 3A 2 — A (3A 3 — A 2 ) = 0. The latter con- 
dition can be understood as a symmetry requirement 
on the potential V by expressing it as a function of the 
fluctuation X\: the most general potential preserving 
rotations, with vanishing V and V', can be written di- 
rectly in terms of the lorentz breaking masses A^: 

V = A tr[A!P t ] 2 - A 2 tr[A 1 P s X 1 P s ] + AgtrLA:^] 2 + 
+2A 4 tr[A 1 P t ]tr[A 1 P s ] + 2AitrLYi P s XiP t ] + • • • (16) 

with P 4 =diag(l/c 2 , 0,0,0), P s = diag(0, 1, 1, 1). In the 
LB phase the last term is absent. Requiring A M = 0, 
we find that the potential is invariant under 



8X1 = ediag(3A 4 c 2 , -A , -A , -A ) 



(17) 



that shifts only the scalar fields ip- and ct_ . This trans- 
formation matches the one encountered for goldstone 
fields in the effective description of [11]. 

However, if we further restrict 3A 4 = — Ao (7=1 
in [11]) we find the "anti-diagonal" Weyl transforma- 
tion 



WeyL 



bXi^tX & X^(l + e)X, (18) 



generated by two opposite rescalings of gi and <? 2 . We 
have thus an explicit non-perturbative symmetry pro- 
tecting [i = 0, that gives in addition A = — 3A 4 . Still 
the graviton mass A 2 can be nonzero and arbitrary. 
One example of a full potential invariant under Weyl_ 
is V = a + aitr^trfA- 1 ] + a 2 tr[A 2 ]tr[A- 2 ] . 
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In the LI phase also there are interesting conse- 
quences: Weyl_ makes all masses proportional, so that 
one cannot get rid of ghosts (A+B = 0) without setting 
all masses to zero. This suggests that the LB phase is 
the only physical one with massive gravitons. 

It is interesting to note that /i = even when the po- 
tential is homogeneous of a generic degree a ^= under 
Weyl_: V(XX) = \ a V{X). Interestingly enough, for 
ansatze of constant curvatures this symmetry leads to 
a constraint on the curvatures as a function of a, with 
possible phenomenological consequences. For instance, 
potentials of degree ±1 lead to vanishing curvature for 
gi : 2- This opens up even the possibility to investigate 
the stability of the LB phase under quantum correc- 
tions effects from matter [10]. 

Phenomenology. First of all it is clear that since the 
lagrangian of our matter C\ is defined only with met- 
ric <7i, no lorentz breaking effects can be observed in 
its propagation and interactions, and the Weak Equiv- 
alence Principle is satisfied [18]. The same happens 
also for matter 2, though the speeds of light in two 
sectors are different. Due to the interaction between 
metrics, the Strong Equivalence Principle is violated 
in both sectors. 

In the weak-field regime, besides the newtonian in- 
teraction oc 1/r we have an additional term that grows 
linearly with r. Moreover, through this term the two 
kinds of matter see each other. For instance, an ordi- 
nary source M will produce potentials in both sectors, 
<Jq = -GM(l/r + n 2 r) and $ 2 = +Kciv 2 GMfi 2 r. No- 
tice that the effect of the linear term is opposite in the 
two sectors. 

Of course if [i = as discussed above this effect is 
absent: both potentials are newtonian, without mix- 
ing between the two sectors. On the other hand if 
/i / the linear term generates a constant acceler- 
ation in direction of the source, attractive or repul- 
sive depending on the sign of /j, 2 [11]. (Amusingly, this 
constant acceleration can explain the Pioneer anomaly, 
a ~ l(T 9 m/s 2 , for fi ~ l(T 21 eV). 



The linearized solution can be trusted only from the 
Schwartschild radius ruv = GM up to distances where 
the linear term drives <J> to be ~ 1; for instance, with 
matter of type 1, rjn = (GM /-t 2 ) -1 ■ In particular for 
the sun r IR ~ 10 12 pcx (^/10~ 21 cV)~ 2 , and clearly the 
linearized approximation works fine for /i ~ 10~ 21 eV; 
on the other hand, taking the galaxy as the source, rm 
drops well below the galactic radius. In such a situation 
a non-perturbative solution is needed [10]. We just 
observe here that since the curvature of the linear term 
vanishes as 1/r, one can hope to match with a well- 
behaved solution at infinity. 

Pulsar binary systems (BPS) set stringent limits on 
the rate of gravitational waves (GWs) emission [13]. In 
our theory the main effect is due to the fact that matter 
in sector 1 emits a combination of graviton mass eigen- 
states, and the massive one is forbidden if the energy 
is too low. The emission rate will thus be modified 
only if m g is higher than the GW energy [19]. For 
BPS a rough limit is m g < 10~ 20 eV. However the BPS 
analysis involves the periastron advance rate [13] that 
is expected to be modified as well in our case, there- 
fore a complete study is needed. Even if m g is below 
such limit, there could be observable effects in GW 
detection: as mentioned, the two graviton states have 
different speeds, therefore it would be possible to ob- 
serve a delay or anticipation of the GWs with respect to 
optical signals coming from the same source. Finally, 
supposing that also type-2 matter can be a source of 
GWs, due to oscillations one could see GW signals in 
our detectors not associated with any otherwise visible 
sources like neutron stars etc. 

We conclude that in this theory the spontaneous 
breaking of lorentz, required to obtain an healthy mas- 
sive gravity theory, turns into an rich phenomenology, 
namely massive, oscillating and lorentz breaking grav- 
itational waves. 
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